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Abstract 

We consider a Dirac one-electron atom placed in a weak, static, uniform magnetic field. We 
show that, to the first order in the strength B of the perturbing field, the only electric multipole 
moment induced by the field in the ground state of the atom is the quadmpole one. Using 
the Sturmian expansion of the generalized Dirac-Coulomb Green function [R. Szmytkowski, 
J. Phys. B 30 (f997) 825; erratum 30 (f997) 2747], we derive a closed-form expression for 
the induced electric quadrupole moment. The result contains the generalized hypergeometric 
function 3F2 of the unit argument. Earlier calculations by other authors, based on the non- 
relativistic model of the atom, predicted in the low-field region the quadratic dependence of 
the induced electric quadrupole moment on B. 

Key words: electric quadrupole moment; magnetic field; Zeeman effect; Dirac one-electron 
atom; Dirac-Coulomb Green function; Sturmian functions 

PACS: 31.15.ap, 32.10.Dk, 32.60.+i 

1 Introduction 

In the mid 1950's, Coulson and Stephen [1] pointed out that the uniform magnetic field should 
induce the electric quadrupole moment (EQM) in the hydrogen-like atom. They used the per- 
turbation theory and found that in the case of the non-relativistic atom, for field strengths B 
corresponding to the complete Paschen-Back effect, the leading term in the expansion of the in- 
duced EQM in powers of the field strength is quadratic in B. For the atomic ground state, their 
result, after being translated into the form conforming both to the definition of the EQM (cf. Sec. 
3 below) and to the notation used in the present work, is 



(2) _ _ 5 0^52 

- 16 Z6 Bi' ^ ' 



where Ze is the nuclear charge, ao is the Bohr radius, a is the Sommerfeld fine-structure constant, 
whereas 

B, = l^^ = f^^6.26T (1.2) 
47r zcOq 
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is the atomic unit of the magnetic induction (pQ is the vacuum permeability and hb is the Bohr mag- 

(2) 

neton). Later, Turbiner [2, Sec. 5] arrived at the same expression for Q20 ^'^^ the Schrodinger one- 
electron atom with Z = 1 and found also explicitly the next non- vanishing term (being Q^q ~ B*) 
in the expansion of Q20 in powers of B. Moreover, using the variational technique, he determined 
the function Q2o{B) for magnetic fields ranging from vanishing to superstrong. A decade ago, 
Potekhin and Turbiner [3] calculated Q2o{B), over a still broader range of B, by two different 
methods, one being the variational approach with a more sophisticated trial function than the one 
used in Ref . [2] , the other being based on the expansion of the perturbed electronic wave function 
in the Landau orbitals. In the low-field limit, the results of Ref. [3] reproduced the quadratic 
dependence of the induced EQM on B predicted in Rcfs. [1,2]. 

In all the aforementioned works, the atom has been described non-relativistically. In the present 
paper, we show that if the relativity is taken into account and the atomic model adopted is the one 
based on the Dirac equation for the electron, it appears that the weak, static, uniform magnetic 
field induces in the ground state of the atom the EQM which, to the lowest order, is linear in 
the perturbing field strength B; no other electric multipole moments are induced in the system 
to the first order in B. Exploiting the Sturmian expansion of the generalized (or reduced) Dirac- 
Coulomb Green function, found by one of us in Ref. [4] and subsequently successfully used in 
analytical calculations of various electromagnetic properties of the Dirac one-electron atom [5-9], 
in Sec. 4 we derive a closed-form expression for the induced EQM, in terms of the generalized 
hypergeometric function 3F2 with the unit argument. 



2 Preliminaries 

It has been already stated in the introduction that the system we shall be concerned with in the 
present work is the Dirac one-electron atom. Its nucleus will be assumed to be point-like, infinitely 
heavy, spinless, and to carry the electric charge Ze. In the absence of external perturbations, the 
atomic ground-state energy level 

^(0)=mcV, (2.1) 

where 

7« = - (aZr, (2.2) 
is two-fold degenerate, the two associated Hamiltonian eigenfunctions, orthonormal in the sense of 

' dV*(o)t(r)vI/(':)(r) = V, (2.3) 

being 

.T.(0)/„^ _ _ , 

r V iQ(°nOf^i^K) 

The radial functions appearing in Eq. (2.4), normalized according to 

'dr{[p(o)(r)p + [0(o)(r)p} = l, (2.5) 



f 

Jo 



are explicitly given by 



^"^^ V«or(27i + i) Uo ; 



(2.6a) 



Q(0)(,) = ./I ( ^) ' e-^^/-, (2.6b) 

^ ^ ^ Y ao r(27i + l)\aQ ) ^ ' 

whereas ^^^(nr), with rir = r/r, are the orthonormal spherical spinors defined as in Ref. [10]. 
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In the presence of a weak, static, uniform magnetic field B = Buz, the level E'^"^ splits into 
two. Their energies, to the first order in B, are given by 

+ 41) (M = ±i), (2.7) 

with 

i?W=sgn(M)^^^A*BB. (2.8) 
The corresponding wave functions, to the same approximation order, are 

*^(r)^*(°)(r) + *W(r) (/. = ±i). (2.9) 

Here, the zeroth-order component ^|f-*(r) is given by Eq. (2.4) (the space quantization axis is 
now chosen to be directed along B). The correction "^'j}\r) solves the inhomogeneous differential 
equation 



s2 

(47reo)r 



7p2 -| ri -| 

-ichcx ■ V + pmc'^ E^°'^ ^'j!\r) = - -eca ■ {B x r) - E^^^ ^'^°\r) (2.10) 



(a and /3 are the standard Dirac matrices), subject to the usual regularity conditions and the 
orthogonality constraint 

/ dVvl/(°)t(^)xi/W(r)=0 {n,n' = ±^). (2.11) 

The integral representation of (r) is 

*(i)(r) = -iecB- / d^r' &°\r,r'){r' xa)'^'^°\r'), (2.12) 
2 Jrs 

where G(°^(r,r') is the generalized Dirac-Coulomb Green function associated with the ground- 
state energy level (2.1). 

3 Analysis of electric multipole moments of the atom in the 
magnetic field 

After these preparatory steps, we set the problem: which electric multipole moments, apart from 
the monopole one, characterize the electronic cloud of the atom in the perturbed state \l/^(r)? If 
Pfj.{r) is the electronic charge density for that state, the spherical components of the L-th order 
electric multipole moment tensor are defined as 



Qlm^ = \I^ZTlL '''^^^MinMr), (3.1) 

where the asterisk denotes the complex conjugation and Y^Minr) is the normalized spherical 
harmonic defined according to the Condon-Shortley phase convention [11, chapter 5]. For the 
atom being in the state ^^(r), the density p^ir) is given by 

-e\l/]',(r)^'„(r) 

P^{r) = • (3-2) 

^3d3r'4(r')*^(r') 

Using Eqs. (2.9) and (2.11), to the first order in the perturbing field, one has 

p^(r)^pW(r) + pW(r), (3.3) 

with 

p(°)(r) = -evl/(0)t(r)vl/(0)(r) (3.4) 
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and 

pW(r) = -e fvI'Wt(^)vi.(o)(r) + «/(")t(r)vl/W(r)l . (3.5) 



Qlm^ ^ Qi°L + Si'L, (3.6) 



Accordingly, it follows that 
where 

Q^M, = -^/^^ *rnOr-^^rVK)*r W (3-7) 
is the multipole moment for the unperturbed atom, whereas 

with 



QiX = -^\/^rhL *r'(^)^'^^^VK)*i.')(r), (3.9) 

is the first-order correction induced by the perturbing magnetic field. To arrive at Eq. (3.8), we 
have exploited the well-known identity 

YLMiUr) = H^Yl^-Mirir). (3.10) 

If the representation (2.12) of ^{/^^(r) is plugged into Eq. (3.9), this gives S^^^ in the form of the 
double integral over M^: 

2iX = \\f^i^"''^ j^f"- j^/^' ^f\ryyiM{nr)G^'\ry)n, • (r' x cx)^f{r'). 

(3.11) 

Using Eq. (2.4) and the explicit representations of the spherical spinors 0,^:\n{nr), it is easy to 
show that the electronic charge density in the unperturbed atom, p'^\r), is spherically symmetric 
and the integral in Eq. (3.7) differs from zero only if L = and M = 0. Hence, for the ground 
state of the unperturbed atom it holds that 

= -e<^io<5Mo, (3.12) 

i.e., all permanent electric multipole moments of the electronic cloud other than the monopole one 

do vanish. 

Next, we proceed to the analysis of the first-order induced multipole moments Q^iXi^- To this 



end, we have to invoke the multipole representation of the generalized Green function, which is 

CO |ft|-l/2 

rr' 



CO \K\-l/-Z 

&^\ry) = ^ E E i 



K=-oo m=— |k|-|-1/2 

9^°}++)ir,r')n^m{nr)nl^{n'^.) -i5^°|_^_) (r, r')0„„(n^)rjl^„(n;) 
i5S_+)(r,r')f^-«™(n.)17L«) 9f\--)ir,r')n_,minr)nl^^in',) 



(3.13) 



If Eqs. (2.4) and (3.13) are inserted into Eq. (3.11), the angular integration in the resulting formula 
may be conveniently carried out with the aid of the identity [10, Eq. (3.1.6)] 



Auk v(k+|)^-A*^ 
riz ■ {rir X (t) f2K^(nr) = i ^^^ _ ^ 0_^^(nr) + i \2k+1\ ^i^+^A'^r) 



V -f2«_i,^(n,) (3.14) 



1^2 _ ^2 



2k- 1 
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and the relation 



2L 



- f d'^nrfliu,{nr)YlM{nr)flK'u.'{nr) 



with 



(3.15) 



In Eq. (3.15), ( ] denotes the Wi gner's 3? coefficient, whereas 

^ ^ \ ma mi, m.c J ° •' ' 



1 

2' 



(3.17) 



and similarly for I' . Exploiting the selection rules embodied in Eq. (3.16) and properties of the 3j 
coefficients, one deduces that the only case when Q^^^f/i does not vanish is the one with L = 2 and 
M = 0. Since Q20/i is real, from Eq. (3.8) one has 

QIm^ = Q^i^L2<5MO, (3.18) 
with Q2o]n being given in the form of the following double radial integral: 

Q« =sgn(M)^(47reo)cS^"dr^°°dr' ( P(o)(r) Q(°)(r) ) r^G^^^r-, r-')r' ( ^|o![^^j ) > 

(3.19) 

where 

Gi°Hr,r')= f Sh!'''!! 1 (3-20) 



is the radial generalized Green function. 



4 Evaluation of the induced electric quadrupole moment 

In the preceding section, we have proved that, to the first order in B, the only electric multipole 
moment induced in the electronic cloud of the atom by the weak, static, uniform magnetic field is 
the quadrupole (L = 2) one, and that its only non- vanishing spherical component is the one with 
M = 0. Now, the Cartesian components of the traceless tensor of the induced quadrupole moment, 

4>=/ dVi(3nr,-r%)p«(r) (z, j G {x, y, 4), (4.1) 

are related to its spherical components through 

o(i) n(i) o<'' \ 

^xxfi ^xyfi ^xzjj, \ 

qW n(i) 

^yxfj, ^HHH ^yzpi 

n(i) qW j 

^zxfi ^zyfj, =iCzzjj, I 

I -in*'' + ./Irn^^' + o''^ 1 



(4.2) 
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Hence, it follows that in the Cartesian basis all the off-diagonal elements of the quadrupole moment 
induced in the system under study do vanish: 



{i j; i,j G {x,y,z}), 



while the diagonal elements are given by 

Q(l) =Q(1) 



2 -^ZZfJ, 



2 ^20^- 



(4.3) 



(4.4) 



Such a structure of the Cartesian representation of the quadrupole moment tensor is characteristic 
for systems possessing the rotational symmetry around the z axis (recall that in our case the z 
axis is directed along the magnetic field). 

It remains to evaluate the double integral in Eq. (3.19). For this purpose, we shall make use of 
the Sturmian expansion of the pertinent generalized radial Dirac-Coulomb Green function for the 
atomic ground state, which is [4] 



G(°)(r,r') 



E l ( Sn}K{r) \ [ (0) „(0) , „(0) / ,x 

. = _Oo MrirK ~ J- \ ^nrK\r) J 



(AC 7^-1). (4.5) 



Here 



/ (l + 7i)(|n^| + 27^)|n^|! 
2ZA/-„,«(iV„^« - K)r(|n^| + 27«) 



(2Zr 



-Zr/ao 



L 



(27.) f'2Zr\ K-^Nn^ (^2^^) f 



2Zr 



:.6a) 



and 



I (l-7i)(|n,| + 27^)|n,|! 

)r{\nr\ + 27„) 



2Zr 
ao 



-Zr/ao 



-(27J 



K- N, 



|nr|-l 



\ ao J \nr\ + 27k 



2Zr 
ao 



(4.6b) 



[with Li"\x) denoting the generalized Lagucrre polynomial [12]; wc define L^^l{p) = 0] arc the 
radial Dirac-Coulomb Sturmian functions associated with the hydrogenic ground-state energy level, 
whereas 

(0) _ i^^rj +7k +iV»^K 



71 



1 



(4.7) 



with 



NurK = ±V(|nr|+7«)^ + (Q:^)^ = ±VKr|2 + 2|n..|7«-hK2 



(4.8) 

being the 'apparent principal quantum number' (notice that it may assume positive as well as 
negative values!); the following sign convention applies to the definition (4.8): the plus sign should 
be chosen for > and the minus one for rir < 0; for rir = one chooses the plus sign if «; < 
and the minus sign if k > 0. Insertion of the separable expansion (4.5) into the right-hand side of 
the formula in Eq. (3.19) leads to the following expression for Q^q^: 



a 



(1) 

20/i 



2 -j^ /"C 

sgn(M)-(47reo)ci3 ^ -(^1 



dr r 



dr' r' 



:(o) 



Mi%Q(°Hr-')4%(r-0+J'^°nr')7^i!kr') 



(4.9) 



The radial integrals in Eq. (4.9) may be taken after one makes use of Eq. (4.7) and of the ex- 
plicit representations of the radial functions P^°\r), Q^'^\r) and the radial Sturmians s'^^J^i''')^ 
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Ti^\{r), given by Eqs. (2.6) and (4.6), respectively. Exploiting the known integral formula [13, Eq. 
(7.414.11)] 

f dx x^e-^L(")(x) = [Re(,) > -1] (4.10) 

and the trivial but extremely useful identity 

7l=7? + 3, (4.11) 

one arrives at 

qW = ccnfu)"''^"" ^ r^(7i+72 + 3) 

^20/. ^8 W ^4 So480(47i + l)r(27i + l)r2(72-7i-2) 



oo 



^ r(|n^| + 72 - 71 - 3)r(|n,| +72-71-2) A^n,2 - 2 

''r.i^^ |n,|!r(|n,|+ 272 + l) 7V„^2 

{\nr\ + 72 - 371 - 3 - 7iA'n.2)(3|n^| + 71 + 372 + 1 + 3jV„^2) 

Krl +72 -71 - 1 + -^nr2 



(4.12) 



The above result may be simplified considerably if in the series X!^=-oo(' ' ' ) ^^^^ collects together 

terms with the same absolute value of the summation index n,. (the Sturmian radial quantum 
number). Proceeding in that way, after much labor, using Eq. (4.8) and again the identity (4.11), 
one finds that 

n(i) = W;.)"'^'^"^ r2(7i+72 + 3) 

^20,. 8 W ^4 Bo 240(471 + l)r(27i + l)r2(72 - 71 - 2) 

r(n,. + 72 - 71 - 3)r(n,. + 72-71-2) 



E 



nr\{nr + 72 — 7i)r(nr + 272 + 1) 



x[(27i - l){nr + 72 - 71 - 3)(n^ + 72 - 7i) + 2(7i + l)(27i + l)(«r + 72 - 7i) 
-6(71 + 1)]. (4.13) 

It is possible to achieve a further simplification. To this end, we express the right-hand side of Eq. 
(4.13) in terms of the hypergeometric functions 2^'i and 3F2 of the unit argument. Since it holds 
that 

r(n + ai)r(n + a2) ^ T{ai)T{a2) f 01,02 
n!r(n + b) T{b) ^ ^\ b 



n}r(r,.^h\ = r(h\ [ b '^J ^^-^^^ 

n=0 

and 



^ r(n + ai)r(n + a2) ^ r(ai)r(Q2) / 01,02,03 A 
^^n\{n + a3)r{n + b) 03^6) V «3 + 1, 6 ' 7 ' 



Eq. (4.13) becomes 



(1) ^ _ . ( s g'eao B r(7i + 72 + 3)r(7i +72+4) 

^20^^ ^4 Bo 1440(71 + 1) (471 + l)r(27i + l)r(272 + l) 

(271 - 1) (72 - 71 - 3) 2F1 f - 7V - ^ 1 



272 + 1 

+ 2(71 + l)(27i + 1) 2F1 ( - - ^ 1 

6(71 + 1) p f 72-71-3,72-71-2,72-71 iA 

72-71 ' 'V 72-71 + 1,272 + 1 'J 



(4.16) 



The two 2^1 functions may be then eliminated with the aid of the Gauss' identity [13, Eq. (9.122.1)] 
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After some simple algebra, wc obtain 



+ 



a^eal B r(27i + 5) 



720(471 + l)r(27i + 1) 
(71 + 72)r(7i + 72 + 3)r(7i +72 + 4) 



r(27i + 5)r(272 + 1) 



2(27^ + 371 + 4) 

72 - 71 - 3, 72 - 71 - 2, 72 - 7i 



72 - 71 + 1,272 + 1 



(4.18) 



Various other equivalent representations of Q2ojj may be derived from Eq. (4.18) with the help of 
recurrence relations obeyed by the 3^2 (1) function. For instance, if one uses repeatedly the relation 



3-f 2 



01,02,03 . 
03 + 1,6 ' 



03 r(6)r(6 - oi - 02) , oi 

H 3-c*2 



oi + 1,02,03 

03 + 1,6 



;i 



oi — 03 r(6 — oi)r(6 - 02) oi — 03 

[Re(6-oi -02) > 0], (4.19) 

and its analogue with oi and 02 interchanged, Eq. (4.18) is casted into the slightly more compact 
expression 

^ a'eal B r(27i + 5) [_ , (71 + l)(7i + 72)r(7i + 72 + 2)r(7i +72 + 3) 
^4 o ...HT^.o- ^ ^+ 7ir(27i + 5)r(272 + l) 



Q 



(1) 

20/n 



X 3-F2 



Bo 1440 r (271) 

72 - 71 - 2, 72 - 71 - 1, 72 - 71 
72 -71 + 1,272 + 1 



(4.20) 



Before concluding, it seems worthwhile to investigate the approximation to Eq. (4.20) for 
aZ <^1. Then one has 

7«^N-^ («eZ\{0}), (4.21) 

and consequently 

r(7« + 7«' + fc) - i\K\ + \K'\ + k-i)\ 













^ i«'i) 









-jiPi\K\ + \K'\ + k) {aZf 

{k,k' e Z\{0}, k G N), 



where 



1 dr(c) 



r(C) dc 

is the digamma function. Using Eq. (4.22) and the recurrence relation 



yields 



V(C + i) = V'(C) + ^ 



(71 + l)(7i + 72)r(7i + 72 + 2)r(7i +72 + 3) ^ ^ lltaZ)^ 
7ir(27i + 5)r(272 + 1) 60^ ' 



Since at the same time, to the second order in aZ, it holds that 

72 - 7i - 2, 72 - 7i - 1, 72 - 7i 



3-f 2 



the sought approximation to QjjOu 



72-71 + 1,272 + 1 
(1) 



1:^1- 



{azy 

40 



20/j 

^(1) / X 23 a^'eon B 

2^°''"^S"(^)240^B^ 



{aZ < 1). 



(4.22) 
(4.23) 

(4.24) 
(4.25) 

(4.26) 
(4.27) 



It is seen from Eqs. (4.27) and (1.1) that for the hydrogen atom {Z = 1) and for the perturbing 
magnetic field comparable with the intra-atomic magnetic field, i.e., for B ~ Bq, the first-order 
quadrupole moment Q^20ix predicted by the relativistic formalism is of the same order of magnitude 

(2) 

as the second-order moment Q20 obtained from the non-relativistic theory. 
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5 Conclusions 



Earlier calculations of the magnetic field-induced electric quadrupolc moment in the ground state 
of the hydrogen-like atom, based on the non-relativistic atomic model, predicted the quadratic 
dependence of that moment on the field strength in the low-field regime. In the present paper, 
we have shown that if the relativity is taken into account and considerations are based on the 
Dirac rather than the Schrodinger or the Pauli equation for the electron, the leading term in the 
expansion of the induced electric quadrupole moment in powers of the field strength appears to be 
linear, not quadratic. The calculations of the actual value of that moment we have carried out in 
Sec. 4 provide a still another example of the usefulness of the Sturmian expansion of the generalized 
Dirac-Coulomb Green function [4] for analytical determination of electromagnetic properties of the 
relativistic one-electron atom. 
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